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O ■ Abstract 

m ■ 

■ New exact solutions for the heat equation with a polynomial non-linearity and for 

the Fisher equation are found. An extended class of non-linear heat equations admitting 

solitary wave solutions is found. The generalization of the Fisher equation is proposed 

43 ■ whose solutions propagate with arbitrary ad hoc fixed velocity. 

Ph' 
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cd ■ 1 Introduction 

a 

The nonlinear reaction-diffusion equations play fundamental role in a great number of 
various models of heat and reaction-diffusion processes, mathematical biology, chemistry, 
genetics and many, many others. Thus, one of the corner stones of mathematical biology 
is the Fisher equation 

Ut - Uxx = u(l - u) (1.1) 

where u = u(x, t) and subscripts denote derivatives w.r.t. the corresponding variable: 
in — dt , u xx — jp • 

Equation (|l.lj) is a particular case of the Kolmogorov-Petrovskii-Piskunov (KPP) equa- 
tion 

u t - u xx = f(u) (1.2) 

where f(u) is a sufficiently smooth function satisfying the relations /(0) = f(l) = 
0, / u (0) = a > 0, f u {u) < a, < u < 1. 

The reaction-diffusion equation with the cubic polynomial nonlinearity 

Ut — u xx = a(u 3 + bu 2 + cu) (1-3) 

where a = ±1, b and c are constants, also has a large application value and includes as 
particular cases the Fitzhugh-Nagumo equation [3] (a = —1,6 = — c— 1,0 < c < 1) which 
is used in population genetics, the Newell- Whitehead 0] (for c = a = — 1,6 = 0) and 
Huxley (for a = b = — l,c = 0) equations. Notice that the Fitzhugh-Nagumo equation 
also belongs to the Kolmogorov-Petrovskii-Piskunov type. 
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A nice property of equations (|l.l|) - ()1.3j) is that they admit plane wave solutions which 
in many cases can be found in explicit form. Existence of such solutions is caused by the 
symmetry w.r.t. translations t — > t + k, x ^ x + r with constant parameters k and r. 
For some special functions f(u) equation admits more extended symmetry groups [5] and, 
as a result, have exact solutions of more general type than plane waves. We notice that 
group analysis of (|1.2|) for f(u) = was carried out by Sophus Lie more than 130 year 
ago 6 . The group classification of systems of nonlinear heat equations was presented in 
papers [Ej. 

The conditional (non-classical) symmetry approach jU], [TQ], enables to construct 
new exact solutions of partial differential equations which cannot be found in the frame- 
work of Lie theory. In application to the equations of type (|1.2|) this approach was 
successfully used for the case of cubic polynomial nonlinearity ()1.3|) only. A systematic 
study of conditional symmetries of equation ()1.3|) was started by Fushchych and Serov 
|12j whereas the most exhaustive analysis of these symmetries was presented by Clarkson 
and Mansfield |13| . In particular, a number of exact solutions for (|1.3|) was found in |12j . 



An effective algorithm for construction of traveling wave solutions together with a 
number of interesting examples was proposed in the recent paper [Tl|. However the 
nonlinear heat equations of the general type <|1-2|) were not analyzed in |14j . 

A goal of our paper is to add the list of known exact solutions of equation (|1.3J1 . Ef- 
fectively we present an infinite number of them. In addition, using the unified algebraic 
method ^1] we select such equations of the type (|1.2|) which admit solitary wave solu- 
tions and construct these solutions explicitly. Finally, we present solutions for the Fisher 
equation and propose such generalization of it which admit the same exact traveling wave 
solution as (1.1), but with any ad hoc given velocity of propagation (for solutions of (1.1) 
this velocity is fixed and equal to -^?). In spirit of Hirota's method 15 to achieve these 
goals we use a special Ansatz which leads to a uniform formulation for all considered 
equations (which, however, is tri-linear). In addition to equations (|1.1|) and ()1.3|) . this 
Ansatz makes it possible to reduce an extended class of equations of the type (|1.2|) . 

In the following section we present a specific Ansatz which will be used to reduce a 
class of nonlinear heat equations. In Section 3 we find an infinite set of new explicit elliptic 
solutions for the heat equation with the cubic and cubic polynomial nonlinearity. 

In Section 4 we describe plane wave solutions for special classes of equations (|1.2|) . In 
Section 5 solitary wave solutions for equations (|1.2j) are found. Finally, in Sections 6 and 
7 we present exact solutions for the Fisher equation and propose a generalization of this 
equation. 



2 The Ansatz and related equations 



We start with the reaction-diffusion equation with a power nonlinearity 



Xu 



n 



A 



2(n + 1) 
(n - l) 2 




where n is a constant, 
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For convenience we choose a special value for the coupling constant A. Scaling u one 
can reduce A to 1 or to -1 for n > 1 and n < 1 respectively. 
For any n ^ 1 we set 

' 7 \ k 9 

, *: = — =-r (2.2) 



2 / n — 1 

and transform (|2.1[) to the uniform equation 

z (z x zt x - z x z xxx - (k - \)z 2 xx ) = z 2 x (z t - (2k + l)z xx ) . (2.3) 

In contrast with (|2.1j) equation (|2.3|) is homogeneous with respect to the dependent 
variable and includes the cubic non-linearities only while (|2.1|) includes u in an arbitrary 
(fixed) power n. We will show that formulation (|2.3f) is very convenient for effective 
reductions. 

We notice that (|2.1j) is not the only nonlinear equation of type (|1.2j) which can be 
reduced to the tri- linear form via Ansatz (|2.2|) . A more general equation (|1.2|) which 
admits this procedure is 



u t - u a 



k (—(k + l)u n + Ain + X 2 u n i 1 + Ag-uV + A 4 ti 2 ~ n ^ (2.4) 



where Ai, . . . , A4 are arbitrary constants. 

Formula (|2.4j) defines special but rather extended class of the nonlinear heat equations, 
which includes all important models enumerated in Introduction and many others. A 
specific presentation of the coupling constants is chosen for convenience. The change 
(|2.2j) transforms (|2.4|) to the following form 



z(z x z x — z x z xxx — \%zz x — A4Z 2 — (k — l)z xx ) 
= z 2 x (z + A12 + X 2 z x - (2k + l)z xx ). 



In contrast with (|2,4[) equation (|2.5|) is homogeneous w.r.t. the dependent variable 
and is much more convenient for searching for exact solutions. 

3 Infinite sets of solutions 

Consider a particular (but important) case of (|2.4|) which corresponds to n = 3, Ai = A2 = 
A 3 = 0: 

u — u xx = — 2u 3 . (3-1) 
The related equation (|2.5|) takes the form 

z(z x 2^xa;a;) — Z x (z &Z XX ). (3-2) 

Using the conditional symmetry approach the following exact solution for (J3.1|) was 
found HB] 

u = 2xds (x 2 + 6t, -^J (3.3) 



where ds(y, k) is the Jacobi elliptic function satisfying 

' drj 
dy 



^l) 2 = k 2 (k 2 -l) + (2k 2 -l) V 2 + V 4 . 



Figure 1: Solution (|3.3|) for equation (|3.1j) 

The plot of this solution for t < 200 created with MATHEMATICA is given in Fig. UJ 

Here we present other elliptic functions solutions for (|3.1[) . effectively an infinite num- 
ber of them. To achieve this goal we exploit conditional symmetry of the potential equa- 
tion (l3~2]) . 

Equation (J3.2j) is compatible with the condition Xz = where 

* = 4-~. (3.4) 
at x ox 

It means |10j that this equation admits conditional symmetry , thus it is reasonable to 
search for its solutions in the form 

z = l p{y), y = x 2 + 6t (3.5) 

where y is the invariant variable for symmetry ()3.4|) , Substituting (|3.5|) into ()3.2|) we come 
to the third order differential equation for (p 

Wyyy = ^VyVyy (3.6) 

Dividing the l.h.s. and r.h.s. of (|3,6j) by (p<p yy and integrating we obtain 

(fyy = C(/9 3 , C = ±2 (3-7) 

where c is the integration constant which can be reduced to 2 (for c > 0) or to -2 (for 
c < 0) by scaling the dependent variable <p. We make such scaling to simplify the following 
formulae. 

In accordance with ()2.2I) . (|3.5|) . any solution ip of (J3.7|) generates a solution for ()3.1|) 
of the following form 

9 TV/}-. 

(3.8) 



2Xfy 



<P 
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An explicit solution of equation (|3.7|) for c = 2 is the Jacobi elliptic function 



r('//i - ds {y,-j=). il r- + 0/ 



(3.9) 



so (|3.8|) can be represented as 



u = U\ 



2xcs (y, ^ j 



(3.10) 



The case c = — 2 leads to the same solution as given in ()3.1U|) . 

The solution (|3.1Ujl is missing in the table of elliptic function solutions present in ^3] • 
The plot of this solution is given by Fig. [2J 




Figure 2: Solution (EHUJ) for equation (jSHJ), < t < 200 

To construct more elliptic function solutions for (|3.1|) we notice that (|3.3|) can be 
rewritten as 

« = yx^(y) (3.11) 

with the same functions ip and y as given in (|3.9|) . Moreover, solution h'J. 1 U\) can be 
obtained from \3.11\) by the change cp — * — . And it is this observation which opens the 
way to construct an infinite number of exact solutions for (|3.1|) . To do this we will exploit 
some properties of the elliptic functions formulated in the following assertions. 

Proposition 1. Let (p = be a solution of equation \3. 7p for c = 2 or c = — 2. 
Then 

(n) 

<P (n+1) = % (3-12) 
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also satisfies this equation for c = 2. 

Proof of this and the following propositions is reduced to a direct verification. We 
notice that equation Q3.7JI is equivalent to the following one 



(n) 



.¥$"J + C " ( 3 - 13 ) 

where C n is the integration constant. Then ip^ n+l ^ of (|3. 12f) satisfies the equation 

(4 n+1 )) 2 = (<4 n+1) ) 4 + C n+1 , C n+1 = -4C n . 

Proposition 2. Xei 6e a solution of equation for C n > 0. T/ien i/tis 

equation is solved also by the following function 



(n) 



(») 



Proposition 3. Lei (/?( n ) be a solution of equation V-i.lty for C n = —B n < 0. Then 
the function 



(n) 



(n) 



satisfies equation \3.1\) for c = —2 and the following relation 



4"> 



+Bl 



(3.14) 



(3.15) 



Using Propositions 1 and 2 and starting with (|3.11|) we obtain infinite sets of solutions 
for equation H3.1|) : 

u n = 2x<p in \ n = 0,l,2,---, (3.16) 



and 



U2k+1 



2 k+1 x 

^(2*+l) 



, fc = 0,l,2,- 



(3.17) 



where {p& k+1 > and <^( n ) are defined by (|3.14|) and the following recurrence relations 

( ra ~i) / i \ 



For n = 0, 1 we obtain from (|3.16j) . (|3.18|) the solutions given by ()3.3|) . (|3.9() . while for 
n = 2, 3, • • • and k = 0, 1, • • • we have 



U2 = 2x 



«3 



2s(«*(y,^)-dn*(y,^)) 



(3.19) 



dn( y ,^)cs( y ,^) (iv^(y,^)-d*(^))' 
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2xdn (y, ^) 




(3.20) 



Formulae (|3.12j) , (|3.19j) , (|3.20|) and the recurrence relations 1)3. 18|) add the list of elliptic 
function solutions for equation (|3.1|) . found in [T5], |13j . Moreover, taking into account 
the transparent invariance of (|3.1|) with respect to displacements of independent variables 
t and x we can write more general solutions changing x — ► x + ki, t — > t + k2 with arbitrary 
constants k\ and k2- 

The plots of solution u\ is given in Fig. |21 



and 




Figure 3: Solution u\ (|3.20j) for equation (|3.1|) 

Propositions 1 and 2 make it possible to construct infinite sets of exact solutions for 
other equations of the type (|2.4|) . 

Consider first equations ()2.4|) for n = 3, A2 = A3 = A4 = 0, i.e., 

u t - u xx = -2 (n 3 + Am) (3.21) 

where without loss of generality we can set Ai = ±1. For Ai = — 1 (j3.21|) is equivalent to 
the Newell- Whitehead [3] equation up to scaling variables t and x. The Ansatze 

u = CxPiOi C = ki cosh(x + ^2) exp(3t), Ai = 1, (3.22) 
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u = rj x (p(r]) , V = ki cos(x + kz) exp(3t), Ai = — 1 (3.23) 

reduces p. 21(1 to the form (|3,7|) . Thus repeating the arguments which follow equation 
(|3.7|) we come to exact solutions for (|3.21|) . The explicit form of these solutions can 
be obtained from Q3.3j) . (|3.16j) - (|3.2U|) Yl& the changes y — ► £, 2x — > £ x for Ai = 1 and 
y — >■ r), 2x -> r] x for Ai = -1. 

Finally, we notice that Proposition 3 makes it possible to construct infinite sets of 
exact solutions for the equations (|2.4|) with n = A4 = — 1, A2 = A3 = 0, i.e., for the 
equations 

iH ~ u xx = 2n 3 (3-24) 

and 

Ut — Uxx = 2 (u 3 + eu) (3.25) 

which differ from l|3.1|) and (|3.21() by the sign of the l.h.s. terms. Indeed, Ansatze (|3.11j) 
and (|3.22|) or 1)3. 23 j) reduce the corresponding equations (|3.24j) and 1)3.25)1 to the following 
equation for 93 

ip" = -2ip 3 (3.26) 

where the double prime denotes the second derivative w.r.t. the corresponding variable 
(i.e., y, £ or rf). 

In accordance with Proposition 3, exact solutions for ()3.26|) have the form (|3.14|) where 
ip( n ^ are solutions of equations (|3.18|) for even n. The related list of exact solutions for 
equation (|3.24|) is given by the following formulae: 



u 



4rrsn (y, ± 

U2 — 



cd (f> 72) ~ dc (v> 72) ~ cn (f> 72) dn ( y ' 72) sn (f ' 72) ' ( ''' 1>] 



2 X 

U2k = ^2*0 



where ip^ 2k ^ are defined by recurrence relations ()3.18|) . 
The plots of solutions 1)3.26)) are given in Figs. El 

Solutions for (|3.25() can be obtained from (|3.27|) by changing y — ► £, x — ► £ x and 
y — > 77, s — ► Tfo for e = 1 and e = — 1 respectively. 



4 Solutions for arbitrary n 

Let us consider equation (|2.4[) with arbitrary n and construct its exact solutions. In this 
section we restrict ourselves to the case A3 = A4 = and use a reduced version of the 
related potential equation ()2.5|) i.e., 

(n - l)z t = (n + 3)z xx - (n - l)(Xiz + \ 2 z x ), (4.1) 

4^2;^ + (n - 3)zl x - (n - l)(Ai^ + \2Z x z xx ). (4.2) 

Any solution of the system (|4.2|) satisfies (|2.5j) with A3 = A4 = 0, the inverse is not 
true. 



S 



2 



Figure 4: Solution u Q (|3.2U|) for equation (|3.1j) 



To solve (|4.1j) . ()4.2I) we introduce the new variable y = Then, dividing (|4.2j) by 
E and using the identity 



z% and using the identity ^f 21 = y x + y 2 we transform this equation to the Riccatti form 



y* + ^y 2 - + Ai) = o. (4.3) 

Differentiating y w.r.t. i and using (j4.1|) and (|4.3j) we obtain the following differential 
consequence 

y = Ay 3 + By 2 + Cy + £> (4.4) 

where 

l(» + S)(»-3)(n+l) B - A 2 (l — — (« 2 - 1 
8 n-1 ' V 16 

C = 1 ((n - l) 2 Al - 2(n + 3)(n - 3)Ai) , £> = ^ T ^A 1 A 2 . 

For arbitrary n the system of equations (|4.3|) . (|4.4|) is compatible but has constant 
solutions only. In three exceptional cases n = ±3 and n = —1 the compatibility conditions 
for (|4.3[) . Q4.4|) are less restrictive in as much as the related coefficient A in (|4.4|) is equal 
to zero. 

Let n / ±3 and -1, then y = c\ = const, and equations ()4.3j) . (|4.4j) . reduce to 
the only condition 

9 2 

Ai = -ciA 2 + {k + l)cf, /c = -. 

n — 1 

The corresponding solution for the system (|4.1j) . (j4.2|) is 



9 



Figure 5: Solution u 2 (|3.20|) for equation (|3.1|) 



and the related exact solution q2.2|) takes the form 



A- ' 



(4.6) 



^1 + C2e -cix-((2fc+l) C 2-A 2 ci)t^ 

Thus we find exact solutions (|4.6|) for the equation 

u t - u xx = -k(k + l)u n + \ 2 ku^~ + ((k + l)c\ - \ 2 ci)ku. (4.7) 
This equation belongs to the Kolmogorov-Petrovski-Piskunov type provided 

A 2 = (fc + l)(ci + l). (4.8) 
The corresponding plane wave solution (|4.6f) propagates with the velocity 

v = k + 1 — kc\. 

In special cases n = ±3 and n = — 1 equations (|4,lj) . ()4.2|) admit more extended classes 
of exact solutions. In particular, for n = 3 we can recover exact solutions for (|1.3[) caused 
by conditional symmetry and classical Lie symmetry as well. We will not study these 
special cases here. 



5 Solitary wave solutions 

Consider now the general equation (|2.3j) with arbitrary parameters X\, X 2 , A3 and A4. 
It seems to be impossible integrate in closed form the related potential equations (|2.5|) . 



10 



Here we search for particular solutions which belong to soliton type and so have good 
perspectives for various applications. 

Let us consider solutions for (|2.3j) of the form z(t, x) = U{£) where £ = fit + x and \i 
is an arbitrary (nonzero) constant. Then we come to the following ordinary differential 
equation for U 

U[U'(pU" - U'" - X 3 U) - X4U 2 - (k - 1)(U") 2 } 

= {U') 2 [{n + X 2 )U' + \ X U - (2k + l)U"} { ° > 



where U' = ^. 

Let us follow |14j and search for solutions for in the form 

U = u Q + vi<p + V2^ 2 H , (5.2) 

where uq, v\, ■ ■ ■ are constants and ip satisfies equation of the following general form 



<p' = e^cq + cup + C2f 2 H (5.3) 

where e = ±1. In order Q5.2|) be compatible with ()5.1|) we have to equate separately the 
terms which include odd and even powers of the square root given by (|5.Hj) . In view of 
this we come to the following system 

U'(jj.UU" - X3U 2 ) = (U' f(v + A 2 ), (5.4) 

U(U'U"' + X4U 2 + {k- 1){U") 2 ) = (U') 2 {(2k + 1)U" - XiU). (5.5) 
Dividing any term in (|5.4|) by fiU 2 U' we come to the Riccatti equation 

Y' - —Y 2 = — 

for Y = jj-, whose general solutions are 

Y = y^tanh(^SS£ + C 



y = y^( tanh (*^F^ + C )) ' if A 2A 3 <0, 



(5.6) 



Sanf^A^ + C'j, if A,A,X). (5.7) 



if A 3 = (5.8) 



A 2 V l l 
/' 



A 2 (£ + C)' 

where C is the integration constant. 

Thus all solutions for (|5.1[) which can be obtained with using the algebraic method 
|14j are exhausted by hyperbolic, triangular and rational ones presented by relations (|5.fij) - 

Solutions (|5.6jl . ()5.7|) and (|5.8j) are compatible with (|5.5j) provided 



H = -A 2 , Ai = -k^, A 4 = (1 - k) f^j , A 2 A 3 / 
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and 

Ai = A 4 = 0, A 3 = 

respectively. Using variables 

t = 7 Tv? > y = 7 X ' a = -^{n - 1), v — t 

\ji — 1J n — 1 A2 

we can rewrite the related equation Q2.4J) as follows: 

u T - u ra = (l + vu 1 ~ n ) (-(n + l)u n + v(n - 3)u + ovHr^j . (5.9) 
The corresponding solutions (|2.2j) . (|5.6() for equation (|5.9|) have the following form 

u = (tanh (ft (y - ^=i) + C?)) ^ , (5-10) 

2 

« = (tanh (ft (y - -^t) + Cj J 1_n (5.11) 

where v < and b = (n — l)y^, 

2 

n-1 



V2 



«=(i/)k=t (tan(6(y--^=t) +C) ) (5.12) 



where ^ > and b = (n — l)w^¥, and 



2 ■ 

2 

in— 1 ' ' - s ' 



if 1/ = 0. 

n-2 



»-2"" J ((n-l)(y--^r + Cl 1 (5.13) 

For ^2 > 1 formula (|5.1U|) presents nice solitary wave solutions which propagate with 
the velocity -7=. In the case n = 2 we come to the bell-shaped solitary wave solution 
which will be discussed in Section 7. 

If ^2 < -1 then JETTO is a singular solution whose physical relevance is doubtful. 
However, in this case equation (|5.9|) admits another solitary wave solutions which are 
given now by relation (|5.11|) . 

We see such solutions exist for the extended class of the nonlinear reaction-diffusion 
equations defined by formula (|5.9|) . 

6 Exact solutions for the Fisher equation 

Let us return to Section 4 and consider in more detail the important case n = 2. Setting 
in 1)4.7}) A2 = 0, ci = — 1 and making the change 

t = 6t, y = V6x (6-1) 

we come to the Fisher equation (jl.lj) for u(r,y): 

U T — Uyy = U(l — U). (6-2) 
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Thus the Fisher equation is a particular case of (|4.7|) and so our solutions (|4.6[) are 
valid for ()6.2|) provided we make the above mentioned changes of variables and set c\ = — 1 
in accordance with ()4.8|) . As a result we recover the well-known Ablowitz-Zeppetella |16| 
solution 

u = 7j. (6.3) 

l + c 2 e^ 6 

This solution can be expressed via hyperbolic functions 

u = u x = j (l - tanh - - c U , c = i In |c 2 |, (6.4) 

u = un = — \\ — coth ( ^_ — — r — c | | (6.5) 
4 V V2^6 12 J J ' ' 

for C2 > and C2 < respectively. 

Taking into account the symmetry of (|6.2|) w.r.t. the discrete transformation u — > 1 — u 
we obtain two more solutions: 113 = 1 — u\ and U4 = 1 — U2- Finally, bearing in mind 
the symmetry of (|6.2ft w.r.t. the space reflection y — > — y we come to four more exact 
solutions by changing y — > — y in Mi, «2i ^3 and U4. 

Thus starting with our general formulae (|4.5j) and (|4.6|) we come to the family of eight 
exact solutions for the Fisher equation. All of them are plane waves propagating with the 
velocity 

To find additional exact solutions we use the Ansatz (compare with IJ3.11J) ) 

u = 3Zy(p(z) (6.6) 

where (p and z are functions to be found. Substituting (|6.6j) into (|6.2j) we come to the 
following reduced equations 

! T = 5ZyV: 2 1 2 (6-7) 

<±6y6yyy 4, yy — ^ y 

and 

ip zz = 3yj 2 . (6.8) 

We see that (p has to satisfy the Weierstrass equation (|6.8|) which we rewrite in the 
following equivalent form 

(p\ = 4<p 3 - C, if = l<p (6.9) 

where C is the integration constant. 

Starting with ()6.6|) and choosing the following exact solutions of (|6.7j) and ()6.9j) : 

/ 1 5 \ 2 . 

* = exp(--^y + -r), <p = — (6.10) 

we come to the Ablowitz-Zeppetella solutions (j6.3j) for the Fisher equation. 

We notice that relations (|6.1U|) present only a very particular solution of ()6.9|) which 
correspond to zero value of the parameter C. In addition, there exist the infinite number 
of other solutions corresponding to non-zero C. The related functions (|6.6|) are: 

u=^z 2 p(z,0,C), z = exp(--^=y+^T + k\ (6.11) 
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where p(z, 0, C) is the Weierstrass function satisfying equation (|6.9f) for 0. 

In order to solutions (|6.6|) be bounded it is sufficient to restrict ourselves to the case 
when ~-^U + | r + k > 0. Such conditions can be satisfied, e.g., for arbitrary positive 
y and negative r and k. The related solutions can be interpreted as ones describing the 
history of the process because the time variable takes arbitrary negative values. The 
graphics of solutions (|6.6|) for some values of the parameter C are given by Figures EHH1 




2 



Figure 6: Solution (j6.11J) with k = 0, C = 10 2 for the Fisher equation ()6.2|) 

Thus the Fisher equation admits the infinite set of exact solutions which include the 
Ablowitz-Zeppetella solutions (|6.3jl and also solutions (|6.11|) enumerated by two parame- 
ters, C / and k. All these solutions are plane waves propagating with the same velocity 

In the following section we consider generalized Fisher equations which admit exact 
solutions with arbitrary propagation velocities. 



7 Generalizations of the Fisher equations 

Let us consider equation (|5.9|) for n = 2, which takes the following form 

u T — u yy = (u + u)(—3u + o\fu — u). (7-1) 

Relation (|7,1|) is a formal generalization of the Fisher equation in as much as in the 
case (7 = (|7.1j) is equivalent to However, for a ^ equation 1)7. 1J) admits soliton 

solutions (j5.1()j) and (or) solutions (|5.11|) - (|5.13j) and so has absolutely another nature then 
(jl.lj) . Nevertheless for small a it would be interesting to treat ()7.1|) as a perturbed version 
of (jnj). 

Consider equation (|7.1|) in more detail. Let v < then scaling dependent and inde- 
pendent variables we can reduce its value to the following one 

v -» v = --, if v < 0. (7.2) 
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Figure 7: Solution ()6.11J1 with k = 0, C = 10 4 for the Fisher equation (|6.2j) 



Setting then 

3 

u = — — u, a = 3e, t = 3r, x = —VSy (7-3) 
we come to the following relation 

u t ~u xx =u^-u + e (^-u^j J. (7.4) 

In the limiting case e — * equation (|7,4|1 reduces to the Fisher equation in the canonical 
formulation 

In accordance with Section 5 equation (|7.1j) admits nice bell-shaped traveling wave 
solution (J5.1U|) which transforms via changes l|7.2|) . (|7.3f) to the following form 

3 

u= -. — -. r (7.5) 

2 cosh 2 (| (x-^+C) 

Consider now equation (|7.1() for va = = and set u = ^. As a result we reduce 
(|7.1jl to the simplest form 

u T — u yy = —u 2 . (7-6) 



The Ansatz 



leads to the following reduced equations 

Zxxx = 0, Z t = KZ XX , K = 5(3 zb a/6). 
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z 



Figure 8: Solution (j6.11|) with k = 0, C = 10 6 for the Fisher equation (J6.2 

Thus we have z = + nt and the related exact solution for (|7.fi[) is 

(3 ± V6)x 2 + 10(12 ± 5V6)t 



3(x 2 + 10(3±^/6)t) 2 



We notice that this solution can be found also using the classical Lie reduction. 

Finally, let us consider one more generalization of the Fisher equation given by relation 
(|4.7j) for n = 2. Using notations (j6.1|) we rewrite it in the following form 

u T — u yy = u(— ci + (ci + l)n2 — u). (7-7) 

For ci = — 1 equation ()7.7|) reduces to the Fisher equation (|6.2j) . 

In accordance with the results presented in Section 4 equation (|7.7I) admits exact 
solutions l|4.6j) which in our notations can be rewritten as 

% = |(l + t a nh(^ y+ ^)r-c)) 2 , 
n 2 = l(l + coth(^ y + Hlfc3) r _ c )) 2 . 

Two more solutions can be obtained by changing y — > — y in (|7.8[l . 

Formulae (|7.8j) present the analogies of solutions (|fi.4|) . (|(i.5j) for equation (|7.7jl . In 
contrast with (|6.4|) . (|6.5() these solutions describe a wave whose propagation velocity is 
equal to 2c ^ 3 Thus changing parameter c\ in ()7.7() we can obtain solutions Q7.8JI with any 
velocity of propagation given ad hoc. In other words we always can take this velocity in 
accordance with experimental data. 

Thus equation (|7.7[) reduces to the Fisher equation if the parameter c\ is equal to — 1. 
Moreover, both equations (|6.2jl and (|7.7|) admit the analogous exact solutions, (|6.4|) . (|6.5|) 
and (|7.8|) . which, however, have different propagation velocities. 
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8 Discussion 



There exist well known regular approaches to search for exact solutions of nonlinear par- 
tial differential equations - the classical Lie approach jS], the conditional (non-classical) 
symmetries method [5], jllj . |1(J| . |12j . generalized conditional symmetries |17j . etc. These 
approaches present effective tools for finding special Ansatze which make it possible to 
reduce the equation of interest and find its particular solutions. 

However, sometimes it seems that the Ansatze by themselves are more fundamental 
than the related symmetries. First, historically, the most famous Ansatze (like the Cole- 
Hopt one for the Burgers equation) was found without a scope of a symmetry approach. 
Secondly, some of Ansatze are effective in rather extended classes of problems character- 
ized by absolutely different symmetries. In addition, in some cases the direct search for 
the Ansatz is a more straight-forward and effective procedure than search for (conditional) 
symmetries. We remind that the conditional symmetry approach presupposes search for 
solutions of nonlinear determining equations which in many cases are not simpler than 
the equation whose symmetries are investigated |18j . 

The present paper is based on using special Ansatze 1)2. 2 Jl , 1)3.11)) , (|6.6|) which have the 
following general form 



where z is an unknown function of independent variables t, x and tp is a function of z. 
The Ansatze (|8.1|l appear to be very effective for the extended class of nonlinear reaction- 
diffusion equations. In particular, they make it possible to find new exact solutions for 
the very well studied heat equations with cubic and quadratic polynomial non-linearities. 
Moreover, such Anaatze can be used to reduce wave equations of another type, e.g., 
hyperbolic equations. We plane to discuss the related results elsewhere. 

We present the extended list of exact solutions for the Fisher equation and for the heat 
equation with the cubic polynomial non-linearity. In fact we present an infinite number of 
different new solutions which are expressed via Weierstrass functions and Jacobi elliptic 
functions. We believe that these solutions will be useful for applications, in particular, 
for analysis of the related boundary value problems. 

We propose a generalization of the Fisher equation which preserves the type of its 
exact solutions, but predicts another propagation velocity. This property differs (|7.7|) from 
numerous other generalizations of the Fisher equation refer, e.g., to and references 
cited therein. 

Finally, we find soliton solutions for a number of nonlinear equations ()2.4j) . To make 
this we use the algebraic method ^1] which however was applied not directly to the 
equation of interest (|2,4[) but to the potential equation (|2.5j) . By this we extend the 
class of non-integrable equations which have soliton solutions to the case of appropriate 
quasi- linear heat equations 1)1.2)1 . 

We stress that all these results were obtained with essential using the Ansatz (|8.1|) . 
It seems to be an intriguing problem to find a regular way for searching such "universal" 
Ansatze. 
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